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Abstract 

We consider a new scenario for supersymmetric decaying dark matter without i?-parity viola- 
tion in theories with goldstini, which arise if super symmetry is broken independently by multiple 
sequestered sectors. The uneaten goldstino naturally has a long lifetime and decays into three- 
body final states including the gravitino, which escapes detection, and two visible particles. The 
goldstini low-energy effective interactions are derived, which can be non-universal and allow the 
dark matter to be leptophilic, in contrast to the case of a single sector supersymmetry breaking. 
In addition, the three-body decay with a missing particle gives a softer spectrum. Consequently, it 
is possible to fit both the e + je~ excess observed by the PAMELA and the e + + e - measurements 
by the Fermi-LAT using universal couplings to all three lepton flavors or 100% branching fraction 
into electrons/positrons, both of which are disfavored in the conventional scenario of dark matter 
decays into two or four visible particles without missing energy. 
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I. INTRODUCTION 



The existence of the dark matter has been firmly established and it constitutes about 23% 
of the total energy density in the universe. The nature of the dark matter is one of the most 
outstanding questions in cosmology and particle physics. Many different types of experi- 
ments are deployed to detect the dark matter and to measure its properties, including direct 
detections from the recoils of the nuclei hit by the dark matter particle, indirect detections 
of the cosmic rays from dark matter annihilations or decays, and collider searches by direct 
production of dark matter particles. Recently, there has been an interesting observation of 
anomalous e + /e~ excess in the energy range of 1 - 100 GeV measured by the PAMELA col- 
laboration which may be interpreted as indirect dark matter signals, coming from dark 
matter annihilations or decays inside the galactic halo. In addition, the e + + e~ spectrum 
measured by the Fermi-LAT experiment between 20 GeV and 1 TeV is harder than that 
inferred from previous experiments j^j, which may also be attributed to the contribution 
from the dark matter. In this paper, we consider a new dark matter candidate which could 
naturally produce the excess of the electron/positron flux observed in these experiments. 

The cosmic positrons are one of the prominent signals for indirect dark matter detections. 
For the most popular dark matter candidate, a weakly interacting massive particle (WIMP), 
electrons and positrons can be produced from annihilations of the WIMPs in the galactic 
halo. However, to account for the PAMELA excess, a large boost factor at the order of 
100 or larger is required to increase the annihilation rate [3|-l6|. In addition, large flux of 
gamma rays will be produced in the dark matter annihilations, which is severely constrained 
by the observed gamma ray spectrum 7H10|. As a result, explaining the PAMELA excess 
by annihilating dark matter has a hard time to satisfy the constraints from the annihilation 
cross section and cosmic gamma ray data. Another possibility to generate the observed 
electron/positron spectrum is that if the dark matter particle is not absolutely stable, but 
decays with a very long lifetime. A small fraction of the dark matter particles has decayed, 
producing electrons and positrons in the decay products ll| . The decaying dark matter has 
an easier time to satisfy the gamma ray constraints, but to explain the PAMELA excess the 
lifetime needs to be of the order of 10 26 — 10 27 seconds UQ, which seems to be an additional 
arbitrary parameter coming from nowhere. 

For the decaying dark matter, it is usually assumed that the symmetry that protects the 
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stability of the dark matter particle is not exact, but violated by some highly suppressed 
interactions. It has been argued that the required lifetime can be obtained from a TeV 
scale particle decaying through dimension-6 operators suppressed by the grand unification 



scale m GUT ~ 2 x 10 16 GeV 



12|. In this paper we consider another possibility that an 



exact symmetry is carried by two sequestered sectors, which interact indirectly only through 
the visible sector (standard model). The lightest particle charged under this symmetry is 
absolutely stable. However, the dark matter is made of the next to the lightest particle 
charged under the same symmetry, which is only approximately stable due to sequestering. 
The dark matter particle decays to the truly stable particle with a long lifetime because 
of the highly suppressed interactions between the two sequestered sectors. The standard 
model (SM) particles produced in the decays can be observed, and could be responsible for 
the anomalies in the cosmic ray experiments. 

We show that such decaying dark matter can arise naturally in the goldstini scenario 
proposed recently 13[. In this scenario, supersymmetry (SUSY) is spontaneously broken in 
more than one sequestered sectors. There is a goldstino associated with the spontaneously 
broken SUSY in each sector. The SUSY in different sectors are connected by supergravity 
and only one combination of the goldstini is eaten and becomes the longitudinal mode of the 
gravitino. The other combinations of the goldstini acquire a mass of twice the gravitino mass 
at the lowest order due to the supergravity effect. Assuming i?-parity is exactly conserved, 
and if the gravitino and an uneaten goldstino are the lightest and the next to the lightest 
supersymmetric particles (LSP and NLSP), respectively, the cosmic electrons and positrons 
can be produced from decays of the goldstino dark matter to the gravitino. If the two SUSY 
breaking sectors only interact indirectly through the visible supersymmetric standard model 
(SSM) sector, the interactions responsible for the goldstino decays are highly suppressed and 
the required lifetime for the observed electron/positron excess can be naturally obtained. 

A distinct feature of this scenario is that the dark matter decays dominantly through 
three-body processes, producing a pair of SM particles and another invisible massive parti- 
cle. Most of the studies of decaying dark matter before assumed that the dark matter particle 
decays through a two-body process to a pair of SM particles or a pair of portals to four SM 
particles without additional missing particles. Some exceptions are in Ref. where the 
three-body decays including a neutrino, as well as from internal bremsstrahlung, are con- 
sidered. The constraint on the anti-proton flux, which shows no excess in the PAMELA 
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experiment 



15j . requires that the decays of dark matter particles dominantly produce lep- 



tons. In the case of two-body decays, the muon and tau final states are preferred [10| and 
the direct decay to the electron and positron pair would give a sharp edge on the energy 
spectrum at half the mass of the dark matter particle, which is not seen by Fermi-LAT. On 
the other hand, the electrons and positrons coming from the three-body decays will have 
a softer and smooth spectrum which may still be consistent with other observations. As 
will be shown, the goldstino couplings to the SM particles, unlike the universal coupling 
of the gravitino, are governed by the fractions of the soft-SUSY breaking masses coming 
from different SUSY breaking sectors for the corresponding superpartners. It is easy for the 
goldstino to have preferential decays to leptons if different superpartners receive different 
soft masses from different SUSY-breaking sectors. 

This paper is organized as follows. In Sec. [TT] we derived the goldstini interactions with 
SM fermions using the method of constrained superfields developed recently by Komargodski 
and Seiberg 16[. From the interactions we can calculate the decay rate of a goldstino to the 
gravitino and a pair of SM fermions. The interactions with other SM fields are collected in 
the Appendix. In Sec. IHII we discuss the model of the decaying goldstino dark matter and the 
parameters which can give rise to the PAMELA signal and satisfy other astrophysical and 
cosmo logical constraints. In Sec. |IV]we perform fits to the electron/positron energy spectra 
observed by PAMELA and Fermi-LAT experiments with the decaying goldstino scenario, 
and identify decay modes and parameters which can be consistent with the observation data. 
We then briefly discuss the collider phenomenology of this scenario. Conclusions are drawn 
in Sec. |V] Throughout this paper, we use "goldstini" when we refer to the goldstino fields 
coming from different SUSY-breaking sectors, and "goldstinos" to represent the plural form 
of the same-species goldstino. 



II. EFFECTIVE INTERACTIONS OF GOLDSTINI 



In this section we derive the low-energy effective interactions of goldstini with two stan 



16|. In 



dard model fermions, by the method of constrained superfields introduced in Ref. 
order to highlight the non-universal nature of goldstini interactions, we begin by reviewing 
the low-energy effective interactions of one goldstino. 



4 



Ref. 16[ considers a system of two chiral superflelds X and Q, 

X = x + V29rj + 6 2 F X , Q = q + V26q + 6 2 F Q , (1) 

interacting through the following Kahler potential K and superpotential W 

K = XX + QQ- ^X 2 X 2 - j^QQXX , W = fX. (2) 

SUSY is spontaneously broken by the F-term vacuum expectation value (VEV) of X, where 
the goldstino resides, while Q is the generic matter field such as the quark or the lepton. 
The c and c terms are included in the Kahler potential to lift the unwanted massless scalars. 
The lagrangian right below the scale A is given by 

C= I d 4 6K+ [ d 2 9W+ [ d 2 9W . (3) 



We are interested in finding the interactions at energies much below the soft SUSY-breaking 
mass scale rn so ft, 

E < m soft < A , (4) 

where the scalar components of X and Q are integrated out. The zero- momentum lagrangian 
is given by 

C = -f + \F X + f\* + \F Q \ 2 -j- 2 \2iF x -rif-^ \qF x + xFq - q V \ 2 , (5) 

which gives rise to the following equations of the motion: 

qF x + xF Q - 7] q = , (6) 
2xF x - V 2 = 0, (7) 

The solutions turn out to be independent of the non-renormalizable couplings c and c. After 
substituting the solutions back into the the chiral superflelds, we obtain the constrained 
superflelds: 

X NL = + V26 V + 6 2 F X , (8) 
zt x 

2 

Qnl = -j=r~ ^F Q + V29q + 9 2 F Q . (9) 
These two superflelds satisfy the constraints: 

X 2 NL = 0, QnlX nl = 0. (10) 



Because the zero-momentum lagrangian in Eq. (j3J) vanishes when evaluated at the solutions 
to the equations of motion, the leading effective interactions involving two goldstinos and 
two q's are obtained from the kinetic term of q: 1 



Ceff = ^(r)q)d»(r)q) + .-- . (11) 



Evidently, the interaction is universal in flavors and only depends on the SUSY-breaking 
scale /. 



Next we consider the "goldstini" scenario 13jj where SUSY is broken independently by 



two sequestered sectors. Matter fields in the SSM may interact with the two SUSY breaking 
sectors only via higher- dimensional operators suppressed by Ai and A2, respectively. The 
Kahler potential and superpotential in this case are 



K = E ( XiXi ~ % x ^ ~ ]^ x ^QQ) +QQ » w = E h Xi 

i=l,2 ^ * ' 1=1,2 



(12) 



The form of the superpotential determines the combination eaten by the gravitino G via the 
super-Higgs mechanism. The eaten goldstino and the uneaten orthogonal combination are 
related to the goldstini of the two sectors by 




(13) 



cos 9 — sin 9 \ I Gl 
sin 9 cos 9 I \ £ 

where Gl is the longitudinal component of the gravitino, and we define 

tanfl=|, feff = y/ ft + ft ■ (14) 

To derive the effective interactions at energies much below m so f t , we follow the same pro- 
cedure as in the single goldstino case to integrate out the scalar components in and Q. 
Furthermore, since we are only interested in the leading-order contribution, the computation 
is greatly simplified if we replace all F-terms by their respective VEV's. In the end the zero 
momentum lagrangian is 

C = E [~fi +\F l + M 2 -^\2i t F l -^\ 2 -^\qF l + x t F Q -q Vl \ 2 ^j , (15) 



1 There is another operator of the form, (qa u q)(d fl r]a 1 ' d^rj) , which is sublcading and only generated at the 
loop- level 
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and solutions to the equations of the motions for the scalars are 



>l2 



2/| ' 



(16) 



/i 2 /A? + /l/Al VA? 



/i . fi 

T2 »7l9 + T2 W 



a; 



mf tan 9 — m\ cot 9 



m~ 



(17) 



where m 2 = / 2 /Af is the contribution from each SUSY-breaking sector to the scalar mass 



of Q and m~ 



m\ + m\. 



It turns out that there are two contributions to four-fermi interactions, in contrast to the 
case of a single goldstino. The first one comes from substituting the solution back into the 
lagrangian in Eq. f JT5|) : 

fhf . i 



c 



(0) 



(q (q 



m 2 ~ 



IV 2 



A? 



(18) 



m + /i a? 

while the second one originates from the scalar kinetic term, d^d^q, which is derivatively 
coupled: 



C 



(i) 

2/ 



-p-d*(G L lt)dHG L q) + -± 



m\i^9-m\ cot 9 



'eff 



eff 



m 2 ~ 



d^ODd^Cq) 



1 / ml tan 9 — m\ cot 9 



j < 



eff 
,(0) 



dJCq)d*{G L q)+h.c. . 



(19) 



Notice that while C 2 j is not derivatively coupled, it only involves the uneaten goldstino and 
not the gravitino. In this sense ( is really a pseudo-goldstino. If we are only interested in the 
decay of the goldstino, Cy obviously does not contribute. For the goldstino annihilations 
or scatterings, however, both contributions are equally important since the two derivatives 
in £2/ P un ou t two factors of order m so f t and we have rrisoft/feff ~ 1/A?- 

Concentrating on Cy , which is relevant for decays of the goldstino, we see that the 
interaction involving only the gravitino is still flavor-universal and insensitive to higher 
dimensional operators in the Kahler potential in Eq. (Tl2|) . while those involving the uneaten 
goldstino are non-universal and do depend on details of the ultraviolet physics. Using this 
effective lagrangian we can compute the decay width of the goldstino into two standard 
model fermions plus the gravitino: 



where N c = 3 for quarks and 1 for leptons, x = /m^, and 



F f (x) = (1 - a; 2 ) (2a; 10 + x 9 - 6x 8 + Qx 7 + 4x 6 + lOQx 5 + Ax A + Qx 3 - Qx 2 + x + 2) 



+60 (x 7 + x 5 ) log 



x 



(21) 



For the benchmark scenario in Ref. 



13| 



rn 



C 



2m dL and F f (l/2) « 0.8 



The effective interaction with two fermions is the most relevant one for the purpose of 
this study. However, for completeness, we present effective interactions of the goldstini with 
other SM particles, such as the gauge bosons and the Higgs fields, in the Appendix. 



III. THE DECAYING GOLDSTINO DARK MATTER 

From the effective interactions derived in Section |TT] we see that if the gravitino is the 
LSP and the goldstino the NLSP, the goldstino only decays through dimensions-8 operators. 



It can be cosmologically stable and play the role of the dark matter [13j. However, since the 
goldstino is not absolutely stable, a small fraction of its relic could have decayed and gives rise 
to interesting astrophysical signals. As mentioned in the Introduction, the recent anomalous 
e + /e~ excess measured by the PAMELA experiment may be interpreted as indirect dark 
matter signals, coming from dark matter annihilations or decays inside the galactic halo. 
The decaying dark matter has an easier time to satisfy the gamma ray constraints but 
requires the lifetime to be of the order of 10 26 -10 27 seconds. This long lifetime could be 
obtained from a TeV scale particle decaying through dimension-6 operators suppressed by 
the grand unification scale mcuT ~ 2 x 10 16 GeV 0, Q. Here we point out that such a 
lifetime can also arise naturally from the goldstino decay. 

Consider that SUSY is broken by two sequestered sectors, SI and S2, independently as 
illustrated in Fig. [TJ through the F-term VEV's , f\ and fi-, respectively. The SSM couples to 
both sectors and receives soft SUSY breaking mass terms through operators suppressed by 
energy scales Ai and A 2 . There is one goldstino from each SUSY breaking sector: rji, i — 1,2. 
One linear combination is eaten and becomes the longitudinal component of the gravitino Gl- 
The other uneaten combination, £, will also acquire a mass due to the super grav ity effect, 



which in the leading order is equal to twice the gravitino mass = 2m 3 / 2 13|. Beyond 



the leading order this relation can be modified 17|. In our discussion we will assume that 



the uneaten goldstino is heavier than the gravitino. Assuming fi ^> / 2 , we have / e g ~ j\ 
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FIG. 1: The super symmetric standard model couples to two sequestered sectors which break SUSY 
independently through the respective F-term VEVs. The SUSY breaking sectors only directly com- 
municate with each other only through the supergravity effect. 



and Gl is mostly rji and ( is mostly 772. From Eq. ( 1T91) we see that the four-fermi coupling 
which governs the goldstino decay into fermions is 

1 f mjtane -m 2 2 cote \ J_ A ™f _ 1^ 1 ™l (??] 

if rh\ and m\ are less hierarchical than ff and /|: m\jm\ 3> To have the goldstino 

mass around the TeV scale as the decaying dark matter, yff\ needs to be ~ 10 11 GeV which 
is the scale for gravity-mediated SUSY breaking. The goldstino decaying lifetime to a single 
SM chiral lepton flavor can be estimated from Eq. fl20]) to be, 




rK4xl026 H^rJ lira?J [wih) [4) l#)J< < 23) 

where m| is the SUSY-breaking mass contribution to the slepton coming from f 2 . We see 
that the necessary lifetime to explain the PAMELA positron excess can be obtained for 
of the order 10 7 GeV, which is suitable for generating a gauge-mediated contribution to the 
soft SUSY-breaking masses in SSM. 

While PAMELA observed an excess in the positron signals, it did not see any anomalous 



excess in the anti-proton signals (15J. This implies that the decays of the jroldstinos should 
mostly produce leptons, with the hadronic channels not exceeding 10% |8j if the positron 
excess is to be explained by the decaying dark matter. From Eq. ff22]) we see that the 
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couplings which govern the goldstino decaying to SM fermions are proportional to fn^ 2 /m~, 
the fraction of the soft SUSY-breaking mass of the corresponding superpartner coming from 
f 2 - Therefore, the decay into quarks can be suppressed if the squarks have a smaller fraction 
of their masses coming from f 2 compared with the sleptons. To satisfy the anti-proton 
constraint it requires 

6 x -f < 0.1 — f =► -4 < 0.13 -f , 24 

\ i J \ i / i t 

where the color factor iV c = 3 is included for the quarks. Such a ratio may result from either 

a smaller m~ 2 or a larger m~. The first possibility could arise if the S2 sector preferentially 

gives SUSY-breaking masses to the sleptons over the squarks. For example, S2 could couple 

to the SSM dominantly through the B — L gauge interaction, then it can give 9 times m\ 

to the sleptons compared to the squarks. In addition, if fi induces rhf through the usual 

gravity mediation, the gravity-mediated contributions to the squark masses are generically 

expected to be much larger then the contributions to the slepton masses due to the running 

contributions from the gluino mass. 2 

In addition to suppressing the decays into quarks, the other decay channels which give 
rise to hadrons such as decays into gauge and Higgs bosons also need to be suppressed. The 
effective couplings between the goldstini and the gauge or Higgs bosons are presented in the 
Appendix. The decay widths into these channels are also controlled by the fractions of the 
soft masses originated from f 2 , and hence could be suppressed in similar ways. This can be 
understood as the uneaten goldstino is mostly composed of r/ 2 in the limit /i > / 2 . In the 
example of B — L mediation from S2, since the SM gauge bosons and Higgs fields do not 
carry B — L charges, decays to these modes can be even further suppressed than decays to 
the quarks. 3 Another possibility is that if SUSY breaking in S2 preserves the i?-symmetry, 
then the couplings to the gauge bosons can be naturally suppressed fl3 |. 

An important difference between the goldstino decaying dark matter and many other 



2 It is also possible that the contributions to some of the scalar mass-squareds from one of the SUSY- 
breaking sector are negative. The goldstino couplings to the leptons can be enhanced if the sleptons are 
light from the cancelation of the mass contributions of the two sectors. 

3 It is worth pointing out that the coupling to two Higgs fields is particularly dangerous since, after the 
Higgs held gets a VEV, there is a corresponding two-body decay into only one Higgs boson plus the 
gravitino. Typically the phase space of two-body decay is larger than that of the three-body decay by 
327T 2 , implying that a strong suppression in the goldstino coupling with two Higgs fields is needed. 
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previously proposed decaying dark matter is that the i?-parity is exact in this case and 
;here is still a missing particle (gravitino) from the three-body decay. Most previous studies 
s[ [ll[ focus on the case of dark matter decaying to two or four SM particles without 
any missing particle (other than neutrinos). Since the energy of the decay products is fixed 
in a two-body decay, if the dark matter decays directly to electrons and positrons, the 
electron/positron spectrum will exhibit a sharp edge at the half the mass of the dark matter 
particle even after propagating through the galaxy. The Fermi-LAT measurement of the 
e + + e~ spectrum does not show any sharp feature below the 1 TeV energy [2J. As a result, 
the decays dominantly to muons and taus are preferred. On the other hand, in our case 
the SM particles from the goldstino decay have a smooth spectrum because they come from 
the three-body decay. The goldstino can decay directly to the electrons/positrons and the 
energy spectrum can still be consistent with the Fermi-LAT result as we will see in the next 
section. 

In the early universe, the goldstinos can be generated from both thermal productions and 
decays of superpartners of SM particles. Requiring the correct relic density for the goldstino 
dark matter turns out to put strong constraints on the reheating temperature Tr. The relic 
from superpartner decays is expected to be dominated by the sleptons since the goldstino 
needs to couple most strongly to the sleptons, and the slepton abundance is less Boltzmann 
suppressed if sleptons are lighter than the other SM superpartners. The decay rate of the 
slepton to a lepton and a goldstino is [3] 

2/ o \ / , — o \ 2 / o \ 

(25) 



mfj tan 0-m| 2 cote 



167T \ U 



cff 




Numerically it turns out to be close to the Hubble scale near the typical freeze-out temper- 
ature of the slepton, 



H(T) 



This implies that a significant fraction of the sleptons has decayed to goldstinos before the 
freeze-out. Above the freeze-out temperature the slepton abundance tracks the thermal 
equilibrium abundance and is exponentially sensitive to the temperature (below the slepton 
mass). As it is well known that the WIMP miracle means that the amount of a WIMP 
particle left at its freeze-out temperature is just about right to account for the dark matter 
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if it survives until today. Therefore, the reheating temperature can not be significantly 
higher than the slepton freeze-out temperature, otherwise there will be too many goldstinos 
coming from slepton decays, which will over-close the universe. This consideration requires 

T» < % ■ (27) 

Goldstinos can also be produced directly in the thermal bath radiation. From Refs. [13I . fl^ 
we see that if the reheating temperature is higher than the goldstino mass, the goldstino 
will be over-produced and the parameters in the range of our interest are clearly ruled out. 
If the reheating temperature is below the goldstino mass, the goldstino production, which is 
proportional to the square of the radiation number density at the high energy tail, will be 
suppressed by the Boltzmann factor exp(— 2777.^/Tr). No over-closure of the universe would 
require 

Tr<^. (28) 

If the slepton mass is not much larger than the goldstino mass, this gives a weaker con- 
straint than the constraint from the slepton decays. On the other hand, there is also a 
lower-bound on the reheating temperature because of the need to produce enough SM su- 
perpartners, whose decays result in the right amount of the goldstino dark matter, as in the 



superWIMP 19[ scenario. This typically requires the reheating temperature to be higher 
than the freeze-out temperature Tp of the lightest observable-sector supersymmetric particle 
(LOSP), which can be a slepton or other superpartners like the neutralino. Combining the 
constraints together, there is only a small window to achieve the right amount of goldstino 
dark matter, 

T F (~ for a weakly interacting LOSp) < T R < Min ^} . (29) 

The fine-tuning required to have the correct relic density seems to be a generic problem for 
models with two goldstini. 

The upper bound of the reheating temperature may be relaxed a little bit if there are 
more than two independent SUSY-breaking sectors. The PAMELA signals may come from 
the decay of a heavy goldstino species to a light goldstino species (remembering that their 
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Given that the reheating temperature is below the electroweak scale in the ran ge of pa rameters we consider 



baryogenesis may require non-thermal production of sphaleron configurations 20|,|2l|, or come from moduli 



decays [22|,|23[ 
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masses can receive corrections to the universal lowest order result). A lifetime similar to 
Eq. fl23l) can be obtained with both \ff[ and ~ 10 9 GeV, while the overall SUSY- 
breaking scale y/ f c s remains at 10 11 GeV due to the presence of additional sectors. Both 
the direct production of the goldstinos and decays from the SM superpartners in the early 
universe are then suppressed by a higher scale (10 9 GeV). In this case, the upper bound of 
the reheating temperature may be raised to around itself. 



IV. ASTROPHYSICAL AND COLLIDER PHENOMENOLOGIES 



A. Indirect detections 



In this subsection we discuss the implications of dark matter indirect detection of the 
scenario considered in Section IIIII As emphasized earlier, this framework differs from the 
conventional decaying dark matter model in that the dark matter in our case dominantly 
decays through the three-body process with a missing gravitino. The resulting lepton energy 
spectrum is softer than that in two-body decays without the missing energy, which allows us 
to fit both PAMELA positron excess and the lack of sharp edge feature in the Fermi-LAT 
e + + e~ measurements at the same time, using final states with electrons. In this work we 
will not be concerned with decays into hadronic final states as well as photons, which are 
assumed to be suppressed. 

We use the Bessel function method of Ref. 24J to calculate the positron flux at the earth, 



due to the decay £ — > Gl£ + £ , with the MED model parameters discussed therein, which 
provide the best fit to the Boron-to-Carbon ratio. For the background fluxes we adopt the 



'model 0" presented by the Fermi-LAT collaboration in Ref. 25|, which are parametrized 
in Ref. {l4j . As for the dark matter halo model, we use the Moore profile in Ref. 26 ]. 
Throughout this section we assume that the dark matter density is p & = 0.3 GeV/cm 3 . 
We perform combined fits to both PAMELA and Fermi-LAT data by varying the decay 
lifetime of the goldstino and the overall normalization of the primary e~ component of 
the background flux, as described in Ref. [14J]. Moreover, since the e + /e~ flux at energies 
below 10 GeV measured at the top of the atmosphere is significantly affected by the solar 
modulation effect, we only use data points above 10 GeV in the PAMELA measurements to 
determine the total \ 2 ■ 



13 




10° io 1 10 Z 10 3 

Energy (GeV) 



FIG. 2: Fits to the positron fraction and the total e + + e~ flux measured by PAMELA and Fermi- 
LAT, respectively, using three-body decay of the dark matter (goldstino) into e + e~ pair together 
with a missing particle (gravitino). Here we assume the mass relation = 2mg . The combined 
X 2 per degree of freedom is 1.6. We also demonstrate the sharp edge in the two-body decay spectrum 
in the Fermi fit. 

In Fig. [2] we show the fit to both the positron fraction of PAMELA and the total e + + e~ 
flux of Fermi-LAT from a 2 TeV goldstino decay into a gravitino and a pair of e + e~ with 
100% branching fraction, assuming the leading order mass relation = 2rriQ L . As can 
be seen, both the rise of positron fraction in the energy regime between 1 and 100 GeV in 
the PAMELA data and the hardening of spectrum at around 400 GeV in the Fermi-LAT 
data can be described by the three-body decay into e + e~ pair plus the missing particle. 
The smooth feature of energy spectrum resulted from the three-body decay is evident in 
the figure. In contrast, we also show in the same figure the fit of a two-body decay of 
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FIG. 3: Fits to the positron fraction and the total e + + e~ flux measured by PAMELA and Fermi- 
LAT, respectively, using three-body decay of the dark matter (goldstino) into l + £~ pair together 
with a missing particle (gravitino). Here we assume the dark matter coupling to all three lepton 
flavors is universal. The combined x 2 P er degree of freedom is 1.5. 

dark matter into e + e~ pair to the Fermi-LAT measurement. The sharp edge at is 
still present even after propagation through the interstellar medium, thus disfavoring this 
particular decay channel as the explanation for the Fermi-LAT measurement [h]] . In Fig. [3] 
we demonstrate that reasonable fits to PAMELA and Fermi-LAT can be obtained if the 
goldstino has universal couplings to all three lepton flavors, which is motivated by the flavor 
changing constraints on slepton masses. 

Some general features of the fit can be understood analytically. The PAMELA data 
show a rising positron fraction above 10 GeV, while the Fermi-LAT measurements suggest a 
hardening feature in the region around 400 GeV. For a three-body decay like ( — > G^t r t~ , 
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if we neglect the mass of the leptons, the maximum possible energy of £ + occurs in the 
configuration when the lepton and the anti-lepton are collinear, q 2 = (pi+ + pe-) 2 = 0, and 
Eg- — > 0. In this case, conservation of momentum gives p^ L = p^ — q and hence 

2 2 

4 m) = < ^ (30) 

in the ( rest frame. The position of the Fermi-LAT plateau set a lower bound on the end 
point of the lepton energy, if we were to explain it with a signal component in the observed 
flux, 

from which we conclude that it will be difficult to fit Fermi-LAT measurement with a dark 



matter mass less than 1 TeV, an observation that has been reached previously |10| . 

The conclusion about the heaviness of the dark matter seems quite robust against different 
choices of background fluxes. For example, the Fermi-LAT collaboration provided two other 
backgrounds, in addition to the "model 0" background adopted in this work, which in fact 
give good fits to their data even in the absence of any exotic sources of e + + e~ flux. These 
backgrounds, model 1 and model 2 in Ref. 25], do not provide good fits to pre- Fermi data 
measured by other experiments. We studied the possibility of using these other background 
to perform the fit to PAMELA and Fermi-LAT. The resulting fits are much worse in the 
case of PAMELA data, since the signal flux that can be accommodated by the Fermi-LAT 
data is too small to explain the rise in the positron fraction. 

The PAMELA anomaly alone can be fitted with a much lighter goldstino, if there are 
other sources which can account for the hardening of the Fermi-LAT spectrum around 400 
GeV. In Fig. [4] it is shown that a 500 GeV decaying goldstino can fit the PAMELA positron 
fraction very well. However, the deficit of the Fermi-LAT spectrum above 100 GeV needs 
to be explained by some other sources, such as those discussed in Refs. 
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B. Collider phenomenology 



The collider phenomenology of the scenario, that the gravitino and goldstino are the LSP 



and NLSP, has been discussed in Ref. 
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29| . Here we only give a brief summary, focusing 
on the parameter region in which we are interested. The SM superpartners produced at the 
collider will cascade-decay down to the LOSP. The LOSP will travel some distance before 



1(3 




IU 1 2 1 

10 10 10 
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FIG. 4: Fit to the positron fraction measured by PAMELA using a 500 Ge V goldstino, and the 
resulting signal flux when added on top of the Fermi-LAT background. In this case, the hardening 
of Fermi- L AT spectrum needs to be explained by some other sources. 

decaying to the goldstino. Therefore, the collider signals depend on which SM superpartner 
is the LOSP. Because the goldstino couples most strongly to the leptons, the most natural 
candidate for the LOSP is the slepton in this case. The decay length of a slepton LOSP is 
estimated to be 

If the long-lived LOSP is a charged slepton, the collider signature is very distinct. Measuring 
the charged track can determine the mass of the long-lived particle. In addition, a large 
fraction of them will decay inside the detector, leaving a displaced kink in the tracking 
detector, which allows a measurement of the LOSP lifetime. On the other hand, if the 
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LOSP is a sneutrino, its decay is invisible and there is no distinct feature other than the 
usual missing energy signals for the SUSY events. 

The LOSP may be other superpartners if their couplings to the goldstino is much sup- 
pressed relative to the slepton couplings to the goldstino. The lifetime of the LOSP in this 
case needs to be longer than that given in Eq. For a colored LOSP (gluino or squark), 
it will hadronize and form i?-hadrons; see Ref. 30[ for a review on the experimental searches. 



A fraction of the i?-hadrons could be stopped in the detector and decay later, resulting in 
distinctive signatures 31] . If the LOSP is a neutralino, it will escape the detector most 
of the time, giving rise to the standard missing energy signals for SUSY. However, a small 
fraction of the neutralinos will decay inside the detector, producing 'J^Z, or h. These de- 
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33j . If the direct 



cays can be discovered if the lifetime is shorter than 10~ 3 -10 -5 second 
LOSP coupling to the goldstino is highly suppressed, the LOSP decay to the goldstino may 
be dominated by 3-body process through the off-shell sleptons. If this indeed happens, it 
provides a nice check that the goldstino couples mostly to the sleptons. 



V. CONCLUSIONS 



In this work we proposed a new scenario for supersymmetric decaying dark matter in 
theories with goldstini, where the uneaten goldstino dominantly decays into gravitino, which 
shows up as missing energy, and two SM particles. In this scenario it is not necessary to 
introduce .R-parity violations, since the goldstino decays through dimension-8 operators and 
naturally has a long lifetime suitable to explain the positron excess observed by the PAMELA 
collaboration. We derive low-energy effective interactions of the goldstini and show that the 
couplings can be non-universal, while the gravitino coupling remains universal as expected. 
The non-universality of the goldstini coupling is crucial for the dark matter to be leptophilic, 
so as to avoid the lack of excess in the anti-proton spectrum measured by PAMELA. To 
obtain the correct goldstino relic density for the dark matter, however, seems to require 
some fine tuning of the reheating temperature in the early universe. 

A distinct feature of this scenario is the three-body decay of the dark matter, which 
results in softer energy spectra for the electrons and positrons, as opposed to a sharp edge 
in the case of the more conventional two-body decay. Consequently, it is possible to fit both 
the positron excess in the PAMELA data and the hardening feature in the e + + e~ flux 
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measured by the Fermi-LAT. We find decays into e + + e~ with 100% branching fraction, 
which is disfavored if the dark matter decays into two or four SM particles, could still provide 
reasonably good fits to PAMELA and Fermi-LAT. In addition, universal coupling of the dark 
matter with all three lepton flavors, which may be favored from other considerations, could 
also fit the data well. 

In this work we have assumed the hadronic decay modes of the dark matter, as well as 
prompt decays into photons, are suppressed in order to satisfy constraints from anti-proton 
and gamma ray measurements. However, it is worth pointing out that most studies on 
these constraints are based on the assumption that the dark matter decays into two-body 
final states, while the decay proceeds through three-body channel with a missing particle 
in our scenario. It would be interesting to re-evaluate these constraints in a more model- 
independent fashion for the case of three-body decays with missing particles. 
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Appendix: Effective Interactions of Goldstini 

Here we present goldstini couplings with the gauge and Higgs bosons without detailed 
derivations, which are beyond the scope of the current work and will be presented elsewhere 



34( | . The goal is to demonstrate that, in the limit j\ ^> /2, these couplings are proportional 
to the fraction of soft masses coming from ji- We use 2-component spinors throughout the 
Appendix. 

Effective interactions of goldstini with two U{1) gauge bosons are derived by using the 



19 



following UV interactions: 



/ ^U + E w x ) WJva + h - c - ' (A - 1} 

J \ i=l,2 % ) 

where the field strength superfield W a is 

W a = -i\ a + L% + a^d,X 6 2 , (A.2) 
Li = €D - iFl . (A.3) 

In the above F(F) = F lw o ilv (F /Jill a :,xu ) . Similar to the case of four-fermi interactions, there 
are two contributions to the two-goldstino and two-photon interactions after integrating out 
the gaugino, arising from the zero momentum lagrangian and the gaugino kinetic term, 
respectively, 



4? = -^r f J en f , VFF V = J- ( + ^ ) fjFFr) , (A.4) 
7 v^Mi + /iA 2 rri\ \ Ax A 2 ' 



27 " 2 f 2 



eff 



. mi tan 9 — m 2 cot 9\ 2 — „ . . 

G L Fa-d[FG L ) (— J rjFa ■ d(Fr)) 



mxtang-^cotg ^ 



(A.5) 



where m; = /j/Aj and = m 1 +m 2 . Again we see that the gravitino effective interaction is 
universal while those involving the goldstino are not. There is also a non-derivative coupling 
for the goldstino. 

There is also a three-point coupling contributing to two-body decays of a goldstino into 
;he gravitino and one massive gauge boson, which is nonetheless suppressed by the D-term 



16l . |35j. It is very small and will not be considered here. 

For couplings with the Higgs bosons we consider the following Kahler potential and 
superpotential: 

K = E ( x ' x - % V E E f 1 - g f x l x ) Hfr* ( A - 6 ) 

i=l ^ ' i=l ot=u ^ * ' 

W = J2 -fiXi + /i (l + j-X^ H u H d , (A.7) 

i=l \ * / 

We will work in the limit Aj ^> A so that Xi = rj 2 /(2f i ) as before. In addition, we define the 
soft masses Bi = dj/j/Aj and m 2 a = g ia f 2 /A 2 such that 

H ~ 5j ~ m so/t and ~ m^ o/4 , (A. 8) 
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and only keep contributions up to 0{m 2 j t / fi). Then the effective interactions relevant for 
goldstino decay into two Higgs bosons are 



C 



(o) 

2ft 



c 



(1) 

2ft 



where 



■G L ( [{m 2 Hu + |/i| 2 ) - 5m 2 d (t)\ - bB\i<\> u 



-u -H- d + h.c. 
1 



(A.9) 



^ { («„ + M 2 ) <Pu - B^ d ) G L j id* {5m 2 u <Pi - 5Biuf> d ) C 

+8, { (tfm^ - SBrity C} ^ ((m^ + |/i| 2 ) 01 - G L 
-\-u -H- d + h.c. , 



if, 



: J>L, 5m 2 =mLtan^-mLcot^, a = u, d 

i 

£ = , 55 = fiitan^ - B 2 cotO . 

i 

It can be seen by the equation of motion for the Higgs fields, 

(m 2 Hu + |/i| 2 ) 4>l - Bfi<f) d = Ucf){ , 
(rn 2 Hd + |/i| 2 ) 4>\ ~ B/i(f>u = D0l , 



(A.10) 



(All) 
(A.12) 



(A.13) 
(A. 14) 



that the above interactions are derivatively coupled, as should be for the gravitino coupling. 
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